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REPLACEMENT  POLICIES  IN  A 


BIVARIATE  FAILURE  MODEL 
by 


Thomas  J.  Tosch 
and 

Paul  T.  Holmes 


ABSTRACT;  A two-component  failure  system  in  which  failure  of 
the  secondary  component  affects  the  residual  life  of  the  primary 
component  is  discussed.  From  two  classes  of  policies  and  two 
criteria,  optimal  policies  for  replacing  the  secondary  component 
are  derived. 


1 . Introduction 

The  bivariate  failure  model  considered  is  a special  case  of 
one  used  by  Tosch  and  Holmes  [2]  (see  also  (31).  Of  the  two 
components,  one  is  considered  essential  to  the  functioning  of  the 
system  and  is  called  the  major  component.  The  other  component  is 
minor  in  that  the  system  can  function  without  it.  The  residual 
lifetime  of  the  major  component  depends  on  whether  the  minor 
component  is  functioning  or  failed.  At  the  failure  of  the  minor 
component  a decision  is  made  whether  to  replace  it  or  not.  This 
decision  is  based  upon  the  particular  replacement  policy  being 
used.  It  is  the  aim  of  this  paper  to  find  policies  that  are 
optimal  in  some  sense. 
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In  section  2,  we  formally  define  the  model  and  give  an  ex- 
ample of  its  applicability.  Further,  the  replacement  policy 
classes  P and  n and  the  criteria  and  r2  for  comparing  policies 
are  defined.  In  sections  3 and  4,  details  of  the  calculations 
needed  to  find  optimal  policies  are  given  for  the  classes  P and 
n respectively.  A numerical  example  is  included.  The  mean  residual 
lifetime  function,  mj(t)  is  discussed  by  examples  in  section  5. 

Some  extensions  of  the  results  of  section  4 are  discussed  in 
section  6. 

2 . Definition  of  the  Model  and  Example 

Label  the  two  components  of  a system  by  A and  B and  let  their 
lifetimes  be  given  by  S and  T respectively.  We  let  A be  the  major 
or  primary  component  of  the  system.  B is  the  minor  or  secondary 
component.  This  means  that  the  system  functions  if  and  only  if 
A functions.  We  assume  that  the  residual  lifetime  of  A depends 
on  whether  B is  functioning  or  not.  If  B never  failed  then  S=X, 
where  X is  a given  random  variable.  At  the  failure  of  B,  however, 

A has  a residual  lifetime  distributed  as  U,  another  given  random 
variable.  It  is  further  assumed  that  X,U  and  T are  mutually  in- 
dependent, positive  random  variables  with  finite  means.  A generali- 
zation of  this  model  which  includes  a change  in  the  residual  life 
of  B at  a failure  of  A is  discussed  in  Tosch  and  Holmes  [2] . 

When  B fails,  it  can  be  replaced  instantaneously  so  that 
there  ii  effect  on  the  residual  lifetime  of  A.  If  B is  not 

replacec  • the  residual  lifetime  of  A is  U.  When  A fails. 


the  system  fails  and  all  failed  components  are  replaced.  The 


system  then  begins  functioning  again. 

Let  be  the  failure  times  of  the  system.  Let 

the  minor  component  have  an  exponentially  distributed  lifetime, 
that  is,  T'exp(6).  At  the  failure  of  A,  if  B is  working,  the 
residual  lifetime  of  B is  exponentially  distributed  with  parameter 
6.  Thus  {Sj^ , S2 » • • • • } is  a renewal  process.  Let  S be  a random 
variable  with  the  same  distribution  as  the  inter-renewal  times 
Sj^,S2-Sj^ , . . . . Under  any  replacement  policy  we  have 

E(S)  < E(X)+E(U)  < ® (2.1) 


The  inter-renewal  periods  will  be  called  cycles. 

As  an  application  of  this  model,  let  A and  B be  two  electrical 
generators  supplying  a hospital.  Suppose  A supplies  3/4  of  the 
electricity  while  B only  supplies  1/4.  When  generator  B fails,  A 
is  powerful  enough  to  supply  all  of  the  power  but  with  added  risk 
of  failure  itself.  B,  however,  is  unable  to  supply  all  of  the 
power  needed,  so  that  A must  be  repaired  whenever  it  fails.  When 
generator  B fails,  we  make  the  decision  either  to  fix  it  immediately 
or  wait  until  A fails  to  repair  them  both. 

The  following  notation  will  be  used: 
c — cost  of  replacing  component  B 
K — cost  of  replacing  component  A 

(2.2) 

(R;P) — the  length  of  the  interval  R using  replacement  policy  P 
C(R;P) — the  total  costs  incurred  during  interval  R using  policy  P. 
R may  be  random  or  fixed.  We  now  introduce  the  criteria  and 


( 
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replacement  policices  to  be  considered. 

A natural  criterion  to  use  is  that  of  minimizing  the 
stationary  average  costs  per  unit  time.  That  is,  minimize 

li„  gCijO^tllP)  , 

where  P comes  from  a particular  class  of  policies. 

Since  {Sj^,S2»...)  is  a renewal  process  and  E(S)  * <*>,  by 
a theorem  due  to  Johns  and  Miller  (c. f . , Ross(l),  p.  52),  we 
have 

, EC(  (0,t]  ;P)  EC((0,S1;P)  . (2.3) 

Ji; 1 — = -EiiSiiiTpT 

This  is  the  basis  for  the  following  criterion. 

Definition  1.  Let  be  the  criterion  which  judges  policy  P 
better  than  policy  Q if  and  only  if 

EC( (0,S] ;P)  < EC( (0,S1 ;Q)  . 

E(  (0,S] ;P)  E((0,S];Q) 

It  is  also  reasonable  to  consider  money  already  spent  as 
gone  and  only  consider  future  costs.  In  particular,  consider 
only  the  costs  from  time  t,  a failure  time  of  B,  until  the  end 
of  the  present  cycle. 

Definition  2.  Let  r2  be  the  criterion  which  judges  policy  P 
better  than  policy  Q if  and  only  if 
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EC( [t;Sl ;P)  < EC( [t;Sl ;Q  , 
E(  [t;Si ;P)  ET(t;S] ;Q) 


given  a failure  of  B at  time  t. 

We  will  now  define  the  policies  to  be  considered. 

Definition  3.  Let  be  the  replacement  policy  by  which  B is 
replaced  anytime  it  fails  before  fixed  time  t,  measureed  from 
the  beginning  of  the  cycle.  It  is  not  replaced  if  a failure 
occurs  after  time  t.  Let  P be  the  class  of  such  policies. 

Each  time  B fails,  a decision  is  made  whether  to  replace  it 
or  not.  Of  course,  once  a decision  is  made  not  to  replace  B, 
there  are  no  more  decisions  to  be  made  in  that  cycle.  The 
criteria  selected  require  the  calculation  of  future  costs  in  the 
cycle.  Therefore  to  judge  a policy  we  have  to  specify  what  actions 
will  be  taken  in  the  future  of  the  cycle  and  not  just  at  a failure 
time  of  B. 

Definition  4.  Let  IT  be  the  class  of  policies  where,  given  a 
failure  of  B at  time  t,  our  decision  whether  to  replace  B or  not 
is  based  on  a comparison  between  the  policies  of  not  replacing 
B and  continuing  to  replace  B until  A fails. 

3 . Optimizing  in  class  P with  criterion 


We  seek  to  find  t > 0 so  that  ( (0 ,Sl ; P^)  ^ minimum. 

E( {0,S] ;P^) 


Let 


d{t)=E{  (0 


E((0,S]  ;P^|x=x)dFjj(x) 


= E((0,S];P. |X=x)dF^(x)+ 

Jq  ^ ^ 


t+ 


E( (0,S) ;P^|X=x)dF^(x) . (3.1) 
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On  the  interval  [0,tl,  E ( (0 ,Sl ; | X-x) »x,  since  A fails  before 
we  would  have  stopped  replacing  B.  On  (t,*),  we  have 


f E(  (0,si  ;P.  |x=x)dF^(x)»[  [ E ( (0  ,S ) ; P.  | X«x  ,T-y ) dP  (x)  6e"'^''dy 

Jt+  ^ JQ  Jt+  ^ ^ 

nt+y 

E(  (0,S1  ;P.  |x»x,T-y)dF^(x)  e '’^dy 

. t+ 

+ [ 1 E(  (0,S1  ;P^  |X=x,T«y)dF^(x)  Be'^'^dy  (3.2) 

JoU+y  ^ 


In  both  of  these  integrals  X > t,  so  that  at  time  t both  components 
are  functioning  and  the  residual  lifetime  of  B is  exponentially 
distributed  with  parameter  B.  We  condition  on  this  lifetime  uo 
that  B would  fail  at  time  t+y.  In  the  first  integral  of  (3.2) 
we  have  the  relation  t < x < t+y  which  means  A fails  sooner  after 
t than  B would  have.  This  implies 


I E( (0,S] ;P^|x=x)dF^(x) 

t+y  • GO  mCO 

xdF  (x)6e"®ydy+  (t+y+E(U))dF  (x)6e”®^dy 

. t+  J 0^ t+y  ^ 


x-t 


Be  ®^dyxdFjj  (x)- 


t+y 

00 

(t+y+E(U) ) 6e"®ydydF^(x) 


ff 

J tJ  0 


f oo 

=1 

• QO 

J t 


■C' 


,-B(x-t)  ,^1  _ ^-6  (x-t) 
6 


xdF„{x)  + | [ (t+E(U)  ) (1-e  ‘")+^-e  (x-t+^)  ] dF^  (> 


tx-t-E  (U)  -x+t-i)  +t+E  (U)  +|}dFj^  (X) 


= (t+E  (U)  +|)  Fx  (t)  - (E  (U)  +|-) 


g-6(x-t)^p  (x) 

X 


Substituting  these  quantities  into  equation  (3.1)  to  obtain  an 
expression  for  d(t),  we  have 


d(t)=E((0 


Let  n(t)=EC((0 


< 


.S];P  )=[ 

^ ■»  0 


xdFj^(x)  + (t+E(U)+|)Fj^(t) 

-(E(U)+|)J  e"^^’‘"^’dFj^(x) 


(3.3) 


EC( (0,SJ ;P^|x-x)dF^(x) 


EC(  (0,S)  ;P^|X»x)dF^(x) +J  EC  ( (0  ,S  1 ; P^.  | X«x)  dF^  (x)  . (3.4) 

On  the  interval  (0,t),  EC ( (0, Si ;P^ I X=x) »c6x+K  since  x is  the 
length  of  the  cycle  and  6x  is  the  expected  number  of  times  that 
B had  to  be  replaced.  The  second  integral  in  (3.4)  is  handled  as 


above  in  the  derivation  of  d(t).  After  simplifying  we  have 


n(t)=c(6t+l)F„(t)+K+cS  x.iF„(x) 


[ x.iF  (x)-c[  ' 
J 0 ^ } t 


e"®(x"t)dF^(x) 


(3.5) 


We  wish  to  minimize  for  t ^ 0.  ^ is  differentiable, 

so  that  the  global  minimum  will  occur  at  the  extremes  or  at  a 
time  t where  the  derivative  is  zero.  The  values  at  the  extremes 
are  given  by 


h(0)  K+c(l-f^(6)) 

d(0) 


(E(U)+|)  (l-f*(B)) 


(3.6) 


where  f is  the  Laplace-Stielt jes  transform  of  F , and 

X X 

n(t)  _ K+c6E(x) 

dTtT  Eon  * 


(3.7) 


The  numerator  of  the  derivative  controls  the  sign.  Let 
r(t)  be  the  numerator  of  the  derivative,  then  it  can  be  shown 
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that : 


r(t)=c6E(U)  lF?{t)  + (Bf  xdF^(x)-Fj^(t)+BE(U)tFj^(t) ) J e"® dF^^Cx)  ) 


(B[^ 


^OD 

-K(F  (t)-(BE(U)+l)  I e'^'^^'^^^dF^Cx)  ]. 

X j ^ X 


(3.8) 


This  expression  is  difficult  to  work  with  in  general,  but 
may  be  evaluated,  at  least  numerically,  for  particular  applica- 
tions. 

Example  1.  Consider  the  application  in  the  beginning  of  the  re- 
port. Let  T'exp(4),  X'exp(2)  and  c=1000,  K=3000,  E(U)».25. 

The  graph  of  , 0 £ t £ 3 is  given  below. 


11,000  - 

10,500 

10,000 

0 12  3 

Figure  1.  A graph  of  n(t)/d(t). 

In  this  example,  ^ is  monotone  decreasing,  so  that  the  optimal 
policy  is  to  always  replace  B when  it  fails. 
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The  difficulty  in  this  derivation  probably  lies  in  the  choice 
of  the  class  of  policies.  We  are  assuming  that  if  it  is  not  optimal 
to  replace  B at  time  t+s^,  then  it  is  not  optimal  to  replace  B had 
it  failed  instead  at  t+S2*  where  S2  > Sj^  > 0.  If  the  failure  rate 
function  of  X,  h (t)*f^(t)/F  (t) , is  like  that  in  Figure  2,  that 
assumption  would  certainly  not  be  justified. 


Figure  2.  Example  of  the  Function  hj^(t). 

Here  it  may  not  be  optimal  to  replace  B if  it  fails  at  time 
t+Sj^,  since  A is  approaching  a very  critical  period  where  it  is 
likely  to  fail  anyway.  At  time  t+S2»  however,  component  A is 
very  reliable  and  therefore  it  may  be  optimal  to  replace  B at 
that  time. 

4 . Optimizing  in  the  class  IT^ 

Throughout  the  rest  of  the  paper,  we  assume  that  a failure 
of  B has  occurred  at  time  t.  All  probabilistic  statements  are 
conditioned  on  that  assumption.  The  class  11  and  criterion 
will  be  considered  first.  At  time  t,  component  B has  failed  and 
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we  must  decide  whether  to  replace  B or  not.  We  won't  replace  B 
if  and  only  if  the  policy  of  not  replacing  B is  "better"  than 
always  replacing  B until  A fails.  "Better"  here  depends  on  the 
criterion  being  used>  in  this  case  P^.  Let  be  policy  of  not 
replacing  B at  time  t,  and  11^  be  the  policy  of  replacing  B until 
A fails.  The  optimal  decision  will  be  to  replace  B at  time  t 
if  and  only  if 


EC((t;S);n^)  EC(It,Sl;nQ) 
E(  (t;S)  TrrmTTlT^)’ 


Under  S-t=u  so  that 

E( lt,S) ;nQ)=E(U) . (4.2) 

The  costs  involved  are  those  of  replacing  B and  A once  each  so 
that 


EC( [t,S) ;nQ)=c+K 


(4.3) 


Under  the  residual  length  of  the  cycle  is  the  residual  length 
of  X after  t.  The  expected  value  of  this  length  is 


E( [t,S] ;n^)=m^(t)= 


Fjj  (x)  dx 
1 


F 


X 


(t) 


(4.4) 


I 


L 


I 
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Finally,  the  expected  costs  incurred  during  this  period  are 
c+K+cBm  (t),  the  last  term  enters  in  since  6m  (t)  is  the  ex- 

X * 

pected  number  of  times  we  will  have  to  replace  B in  the  interval 
(t,S].  Therefore 

EC(  lt,S]  ,nj^)=c+K+cBmj^(t)  (4.5) 


Combining  these  into  (4.1)  and  doing  the  algebra,  we  obtain 
Theorem  6 : In  class  H with  criterion  it  is  optimal  to  replace 

a failed  component  B at  time  t if  and  only  if 


(c+K-c6E(U) )m^(t)  > (c+K)E(U). 


Since  ro^(t)  is  not  monotone  in  general,  these  policies  will  not 
have  the  cutoff  point  between  replacing  and  not  replacing  as  did 
the  policies  presented  in  the  previous  section.  The  following 
corollary  is  immediate. 

Corollary  7;  If  c6E(U)  > c+K  then  it  is  optimal  never  to  replace 
component  B. 

G ( U ) K 

The  condition  c6E(U)  > c+K  is  equivalent  to  H '* 

Heuristically  this  says  that  E(U),  the  expected  lifetime  of  A 
without  B,  is  so  much  larger  than  E(T),  the  lifetime  of  B,  that 
it  is  not  worth  the  ris)c  of  c dollars  to  even  have  a component  B 
in  the  system.  Under  normal  circumstances  then,  this  condition 
will  not  be  satisfied. 
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If  X'exp(a),  then  (t) which  implies 
Corollary  8;  If  X'exp(a),  then  in  class  n with  criterion  r2> 
it  is  optimal  to 

i)  never  replace  B,  if  (c+K-BE(U))  £ a(c+K)E(U), 
ii)  always  replace  B,  if  (c+K-6E(U))  > a(c+K)E(U). 

This  result  is  very  intuitive  since  at  any  failure  of  B, 
the  future  looks  exactly  the  same  as  at  any  other  failure  of  B. 
This  implies  that  the  same  decision  should  be  made  at  each  failure 
of  B. 

In  example  5,  we  have  (c+K-c6E (U) ) ^ = (3000) j = 1500,  while 
(K+c)  E (U)  = (4000)  (.  25)  =1000.  So  that  in  11  under  criterion  r^, 
it  is  optimal  to  always  replace  B when  it  fails. 

Basing  the  criterion  only  on  future  costs  has  the  drawback 
that  K,  the  cost  of  replacing  component  A,  is  distributed  only  on 
the  future  time  in  the  cycle.  For  this  reason,  we  will  consider 
the  same  class  II  but  with  criterion  Let  t)e  as  before. 

Under  it  is  optimal  to  replace  B at  time  t if  and  only  if 

EC(  (0,S[;nj^)  ^ EC(  (0,S1  ;nQ) 

~E(  (0,S]  E(  (0,S]  ;nQ 

As  before  it  is  easy  to  obtain  expressions  for  these  quantities. 

We  have 


E( (0,S) ;nQ)=t+E(U) 


(4.7) 


E{  (0,S]  ;n^)=t+mj^(t)  . 


(4.8) 
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These  follow  since  the  cycle  is  already  t units  old  when  we  make 
our  decision.  The  expected  number  of  times  we  have  replaced  B 
in  (0,t)  is  6t  at  a cost  of  cPt  so  that  the  other  quantities 
can  be  written  as: 


EC(  (0,S]  ;nQ)=cSt+c+K,  (4.9) 

EC(  (0,S1  .-II  >=c3t+c+K+cfim  (t)  . (4.10) 

1 ^ 

Combining  these  into  (4.6),  our  decision  is  to  replace  if 

c6t+c+K+cBm^  ( t ) cBt+c+K  , ^ j 

t+nTTtl  t+E  (U) 

from  which  we  get 

Theorem  9 1 In  class  n with  criterion  it  is  optimal  to  re- 

place component  B at  time  t if  and  only  if 

(K+c-cPE(U) )m^(t)  > cBE(U) t+(K+c)E(U) . 

Corollary  10;  If  K+c  < c6E(U)  then  it  is  optimal  to  never 
replace  B. 

The  proof  is  trivial.  The  remarks  on  the  condition  following 
corollary  7 are  valid  here  also. 

Returning  to  example  5,  (c+K-cBE  (U)  ) m^^  (t)  =1500  for  all 
t ^ 0.  (K+c) E (U) =1000  and  cBE(U)=1000.  Therefore  it  is  optimal 
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under  to  replace  B if  and  only  if  1500  > lOOOt+1000,  i.e., 

if  t is  in  the  interval  [0,'j). 

There  is  a large  similarity  in  the  conditions  needed  for 
replacement  under  and  r2*  In  fact  , we  have 
Corollary  11;  If  the  optimal  decision  under  is  to  replace 
B then  it  is  also  the  optimal  decision  under  r2* 

The  converse  of  the  corollary  is  clearly  not  true  as  example 
5 illustrates. 

With  criterion  then,  it  is  optimal  to  replace  B when 
the  mean  residual  lifetime  of  A is  greater  than  the  given  linear 
function.  The  next  section  includes  graphs  of  mean  residual 
lifetime  functions  for  various  distributions  to  illustrate  the 


possible  policies. 

5.  The  Function  mj.{t) 

The  following  figures  illustrate  the  various  forms  of  mj.(t). 


15 

In  this  case  there  will  be  a unique  time  t*,  where  we  will  re- 
place B if  it  fails  before  t*,  but  not  afterwards. 


Figure  4.  mj.(t)  for  a Weibull  Distribution  with  6<  1 


In  addition  to  the  previous  policy,  in  this  case,  there  could 
be  that  if  B fails  in  (tj^,t2)  we  will  replace  it; 

otherwise  we  will  not. 
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Fj (t)=. 5e“^^+. 5e'^  . 

The  same  possible  optimal  policies  exist  here  as  in  the  previous 
example. 


Fj.  (t)=0.7e~^  +0.3e~’^ 

Depending  on  the  given  constants,  there  are  a variety  of  optimal 
policies  that  could  occur  here.  One  possibility  is  that  there 
exists  so  that  we  replace  on  the  intervals  [0,tj^), 

(t2»tj)  and  do  not  otherwise. 

6.  Extensions 

The  costs  that  have  been  considered  are  very  simple.  What 
has  been  done  is  easily  extended  to  more  complicated  costs.  The 


I 


I 
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extensions  will  be  applied  to  optimizing  in  the  class  n under 
criterion 

Suppose  that  the  components  are  not  replaced  but  are  repaired. 

Let  K be  a random  variable  which  gives  the  repair  cost  of  component 
A,  and  A be  the  random  repair  cost  of  B.  If  we  let  K=E(K)  and 
c=E(A),  then  the  analysis  is  exactly  as  before  and  the  optimal 
policy  is  given  by  Theorem  9. 

Suppose  instead  that  the  system  has  to  be  stopped  in  order 
to  make  any  replacements,  and  that  a penalty  cost  of  6 is  charged 
for  each  unit  of  time  that  system  is  not  functioning.  Let 
be  the  random  length  of  time  needed  to  replace  A and  H2,  the 
time  needed  to  replace  B.  Assume  that  the  components  cannot  be 
replaced  simultaneously.  Let  c'=c+yE(H2)/  K'=k+YE{Hj^)  . The 
analysis  goes  as  before  and  the  optimal  policy  is  given  by 
Theorem  9 with  K',c'  substituted  for  K,c. 

If  the  components  can  be  replaced  simultaneously  then  there 
is  a change,  since  if  we  did  not  replace  B when  it  failed,  A 
and  B could  be  replaced  at  the  same  time.  In  this  case,  it  is 
optimal  to  replace  B at  time  t if  and  only  if 

c' 3t+c'+K'+c' 6m„(t)  c' et+c+K+YE(max(H,  ,H-,)  ) ,,  ,, 

X 1 / V o . 1 J 

t+m^(t)  t+E(U) 

The  last  extension  to  be  discussed  deals  with  the  random 


variable  U.  It  has  been  assumed  that  U and  X are  independent. 
This  is  often  not  even  a good  approximation.  Suppose  U is 
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dependent  on  the  residual  lifetime  of  X,  in  that,  E (U) =em^ ( t) +6 , 
e 0,6  ^ 0.  In  this  case,  it  is  optimal  to  replace  B at  time 
t if  and  only  if 


cBt+c+k+cBm^^  (t) 
t+nTTtl 


< cBt+c+K 
t+em^ (t) +6 


(6.2) 


This  reduces  to  the  condition 


cBem^ (t) +m  (t) [cBEt+ (K+c) G- (K+c) +cB6 ] +cB6t+ (K+c) 6 < 0.  (6.3) 

X X 


So  in  this  case,  there  is  a quadratic  relationship  between  t and 
m^(t).  Note  e=0  corresponds  to  the  case  in  section  4. 

Other  extensions  can  be  made  with  little  change  in  the 
form  of  the  optimal  policy.  If  we  do  not  restrict  T to  being 
exponential  then  we  must  include  repairing  B when  A fails  into 
our  model  to  retain  the  renewal  process. 
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